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Introduction tc 
Interplanetary 


■ The interplanetary design problem is composed of both discrete and 

real-valued decision parameters: 

- Choice of destination(s), number of planetary flybys, identities of flyby 
planets 

- Launch date, flight time(s), epochs of maneuvers, control history, flyby 
altitudes, etc. 

■ For example, for a main-belt asteroid mission, the designer must 

choose: 

- The optimal asteroid from a set of scientifically interesting bodies 
provided by the customer 

- Whether or not to perform planetary flybys on the way to the main belt 
and, if so, at which planets 

- Optimal trajectory from the Earth to the chosen asteroid by way of the 
chosen flyby planets 
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Traditional Met 
Flyby Trajector 


■ Several methods of picking the destination and flyby sequence: 

- Grid search over all possible choices of destinations, flyby sequence, 
launch date, etc. (very expensive and often impractical) 

- Intuition-guided manual design of the trajectory (even more expensive, 
can miss non-intuitive solutions) 

■ Several methods of designing the trajectory: 

- Local optimization from an initial guess provided by a chemical 
mission design (but sometimes the optimal chemical trajectory does 
not resemble the optimal low-thrust trajectory) 

- Local optimization from an initial guess provided by a low-fidelity 
approximation to the low-thrust model, i.e. shaped-based methods 
(but sometimes the shape-based method cannot accurately 
approximate the true trajectory) 
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Automated M 
Optimal Cont 


■ Break the mission design problem into two stages, or “loops” 

- “outer-loop” picks sets of destinations, planetary flybys, sizes the 
power system, can pick propulsion system - a discrete optimization 
problem 

- “inner-loop” finds the optimal trajectory for a given candidate outer- 
loop solution - a real-valued optimization problem 

- For the outer-loop to work, the inner-loop must function autonomously 
(i.e. no human interaction) 
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Multi-Object 


■ The customer (scientist or project manager) most often does not want 
just one point solution to the mission design problem 

■ Instead, an exploration of a multi-objective trade space is required 

■ For a typical main-belt asteroid mission the customer might wish to see 
the trade-space of: 

- Launch date vs 

- Flight time vs 

- Deliverable mass 

- While varying the destination asteroid, planetary flybys, launch year, 
etc. 

■ To address this question we use a multi-objective discrete outer-loop 
which defines many single objective real-valued inner-loop problems 
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Outer-Loop T 



■ The outer-loop finds the non-dominated trade surface between any set 
of objective functions chosen by the user 

■ Non-dominated surface means “no point on the surface is superior to 
any other point on the surface in all of the objective functions” 

■ The outer-loop solver may choose from a menu of options for each 
decision variable 

■ The choices made by the outer-loop solver are used to define trajectory 
optimization problems to be solved by the inner-loop 
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I 


Anatomy of 


• Break mission into a set of “journeys,” 
each of which in turn is broken into 
“phases” 


The endpoints of a journey are chosen in 
the problem assumptions 

The endpoints of a phase (i.e. a flyby 
target) may be chosen by the user or an 
Outer-Loop solver 

c 
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Outer-Loop 


Launch Year 
Code Year 

0 2020 

1 2021 

2 2022 

3 2023 

4 2024 

6 2025 


Flight Time Upper 
Bound 


Code 

0 

1 

2 

3 

4 


# Years 

5 

6 

7 

8 
9 


First Asteroid 

ode 

Body 

0 

Ceres 

1 

Pallas 

2 

Juno 

3 

Vesta 

4 

Astraea 

5 

Hebe 


Second Asteroid 
Code Body 

0 Ceres 

1 Pallas 

2 Juno 

3 Vesta 

4 Astraea 

5 Hebe 


First Journey First Flyby 

Code 

Body 

0 

Earth 

1 

Mars 

2 

Jupiter 

3 

No flyby 

4 

No flyby 

5 

No flyby 


7 2026 

8 2027 

9 2028 

10 2029 


Second Journey Flyby 

Code 

Body 

0 

Earth 

1 

Mars 

2 

Jupiter 

3 

No flyby 

4 

No flyby 

5 

No flyby 


5 10 

7 11 

8 12 


6 Iris 

7 Flora 
(475 

^ choices) 


6 

7 


Iris 

Flora 

(475 

choices) 


First Journey Second Flyby 

Code 

Body 

0 

Earth 

1 

Mars 

2 

Jupiter 

3 

No flyby 

4 

No flyby 

5 

No flyby 




Sample Mission 




Flight Time Upper 

Asteroid 1 

Potential Planetary 

Asteroid 2 

Potential 


Bound 


Flyby 1 


Planetary Flyby 2 

Code 

4 

0 

1 

1 

1 

Translation 

By 

Ceres 

Mars 

Pallas 

none 
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Multi-Objecth 


■ The outer-loop optimization problem is solved using a discrete multi- 
objective solver, in this case Non-Dominated Sorting Genetic Algorithm 
II (NSGA-II) 

■ NSGA-II finds the non-dominated front, surface, or hyper-surface 
between any number of objectives chosen by the user 


Initial generation 



Population evolves via 
genetic operators 
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Inner-Loop M 


zation 


■ The inner-loop solves a real-valued trajectory optimization problem 
which is defined by each candidate solution to the outer-loop problem 

■ The inner-loop must function autonomously because the problems are 
generated in real time and there is no opportunity for human 
intervention 

■ The outer-loop is only as good as the solutions to the inner-loop 
problem, so the inner-loop must be robust 

■ A given run of the outer-loop may require hundreds or even thousands 
of runs of the inner-loop, so the inner-loop must be fast 

■ If the individual inner-loop runs are independent then many of them can 
be run in parallel 
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Multiple Gravity 
(MGALT) via the 





Low-Thrust 
agan Transcription 


O patch point 
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Power, Propu 
Modeling 


Medium-fidelity mission design requires accurate hardware modeling 
Launch vehicles are modeled using a polynomial fit 

171 delivered = (l — a Lv) ( a LF^I + b LV C 3 + C LV C f + d LV C f + e LV C 3 + f LV ^j 
where o LV is launch vehicle margin and C 3 is hyperbolic excess velocity 


■ Thrusters are modeled using either a polynomial fit to published thrust and mass flow rate data 

m = a F P 4 + b F P 3 + c F P 2 + d F P + e F 
T = a T P 4 + b T P 3 + c T P 2 + d T P + e T 


or, when detailed performance data is unavailable 

2 riP 
T =7-^ 

*sp do 

■ Power is modeled by a standard polynomial model 



V +Z1 + Z2 
7o + r ' r 2 

1 + y 3 r + y 4 r : 


(1 - zy 


where P 0 is the power at beginning of life at 1 All and r is the solar array degradation constant 


■ Ephemeris data for solar system bodies is provided via the SPICE toolkit 
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Inner-Loop S 
Nonlinear Pr< 


Minimize f(x) 

Subject to: 
x lb <x< x ub 
c(X) < 0 

Ax < 0 

where: 

x ib , x ub are lower and upper bounds on the decision variables 
c(x) is a vector of nonlinear constraints 
Ax is a vector of linear constraints 

■ There are several third party solvers that do this (SNOPT, IPOPT, fmincon, 
vf 1 3AD) 

■ But all of these methods require an initial guess... 


NAVIGATION & MISSION DESIGN BRANCH, CODE 595 

NASA GSFC 


15 * 


Inner-Loop 

Monotonic 



Leary, 2000 

Vasile, Minisci, and Locatelli, 2009 
Yam, di Lorenzo, and Izzo, 2011 
Englander (dissertation), 2013 
Casioli et al., 2013 
Englander and Englander, 201 4 


Improved from standard MBH by: 

1 . “Feasible point finder” aggregate penalty method 

2. Non-uniform (Pareto) perturbation step 

3. “Time-hop” operator (Casioli et al.) 
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Example: Mai 


Mission Objective 

Visit two main-belt asteroids with diameter greater than 50 km 
(475 bodies meet this filter) 

Launch Vehicle 

Atlas V 401 

Power System 


Array power at 1 All 

15 kW 

Cell performance model 

1/r 2 

Spacecraft bus power 

800 W 

Power margin 

15% 

Propulsion System 

Thruster 

NEXT (throttle table 11, high-lsp mode) 

Number of thrusters 

1 

Duty cycle 

90% 

Propellant tank 

unconstrained 

Mission Sequence 

up to two planetary flybys are permitted before the first asteroid and up to 


one between the first and second asteroids 

Inner-Loop Objective Function 

Maximize delivered mass to second asteroid 

Outer-Loop Objective Functions 

Delivered mass to second asteroid 


Launch year 


Flight time 
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Main-Belt T 
Outer-Loop 


Launch Year 

Code 

Year 

0 

2020 

1 

2021 

2 

2022 

3 

2023 

4 

2024 

6 

2025 

7 

2026 

8 

2027 

9 

2028 

10 

2029 


Second Asteroid 

Code 

Body 

0 

Ceres 

1 

Pallas 

2 

Juno 

3 

Vesta 

4 

Astraea 

5 

Hebe 

6 

Iris 

7 

Flora 


(475 

choices) 


First Asteroid 

Code 

Body 

0 

Ceres 

1 

Pallas 

2 

Juno 

3 

Vesta 

4 

Astraea 

5 

Hebe 

6 

Iris 

7 

Flora 


(475 

choices) 


Flight Time Upper 
Bound 

Code 

# Years 

0 

5 

1 

6 

2 

7 

3 

8 

4 

9 

5 

10 

7 

11 

8 

12 


First Journey First Flyby 

Code 

Body 

0 

Earth 

1 

Mars 

2 

Jupiter 

3 

No flyby 

4 

No flyby 

5 

No flyby 


First Journey Second Flyby 

Code 

Body 

0 

Earth 

1 

Mars 

2 

Jupiter 

3 

No flyby 

4 

No flyby 

5 

No flyby 


Second Journey Flyby 

Code 

Body 

0 

Earth 

1 

Mars 

2 

Jupiter 

3 

No flyby 

4 

No flyby 

5 

No flyby 


1.16x10 9 
possible 
combinations, 
4.82x1 0 9 with 
duplicates 
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Delivered mass 


Main-Belt Tw< 
First Generat 



CODE 595 
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Delivered mass 


Main-Belt Tw< 
Final Generat 
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Main-Belt Tw< 
Example Traj 


unpowered flyby 

Mars 

6/9/2025 

= 5.593 km/s 
DEC = 20. 9 C 
altitude = 300 km 
m = 2680 kg 


LT rendezvous 
20Massalia 
7/22/2033 
m = 2055 kg 


LT rendezvous 
554Peraga(1905PS) 
6/18/2030 
m = 2191 kg 



unpowered flyby 
Earth 

9/23/2024 
Vqq — 5.412 km/s 
DEC = -5.2° 
altitude = 300 km 
m = 2779 kg 

launch 

Earth 

7/23/2023 

C3 = 1.311 km 2 /s 2 

DLA = -9.9° 

m = 2951 kg 

departure 

554Peraga(1905PS) 
6/18/2031 
m = 2191 kg 


A 10-year mission launching in 2023 delivers 2055 kg to Peraga and Massalia 
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Conclusions 


■ The low-thrust interplanetary mission design problem may be posed as a multi- 
objective hybrid optimal control problem 

■ The combination of a multi-objective discrete NSGA-II outer-loop with a 
MBH+NLP inner-loop is a very powerful way to explore a mission trade space in 
an efficient, automated manner 

■ The algorithm described here has revolutionized the low-thrust interplanetary 
mission design process at NASA Goddard Space Flight Center 

We can now study multiple mission design cases simultaneously, limited only 
by available computing power 

Mission design engineers can now spend more time with the customer and 
with spacecraft hardware engineers so that we can fully understand the 
scientific and engineering context of our work 

Good mission ideas are much less likely to be rejected due to lack of time to 
work on mission design, and bad ideas are much more likely to be rejected 
before they consume too many resources 
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EMTG is available open-source at 
https://sourceforae.net/proiects/emta/ 
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